Phase Transition in a Two-level-cavity System in an Ohmic Environment 
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We propose that in the presence of an Ohmic, de-phasing type environment, a two-level-cavity 
system undergoes a quantum phase transition from a state with damped Rabi oscillation to a state 
without. We present the phase diagram and make predictions for pump and probe experiment. 
Such a strong coupling effect of the environment is beyond the reach of conventional perturbative 
treatment. 



A combined two-level system and an cavity electromag- 
netic mode has been proposed as an important realiza- 
tion of qubit in quantum information science 0, Q, [H, Ell- 
in quantum optics the coherent oscillation, or Rabi os- 
cillation, of such a combined system is the key to phe- 
nomena such as lasing and electromagnetically induced 
transparency 0, In both cases the environment acts 
to destroy the desired quantum coherence @, S| ■ As a re- 
sult, understanding the coupling to the environment and 
controlling it has become one of the most pressing task 
in quantum control. 

In conventional wisdom, the environment broadens 
quantum levels and damps the Rabi oscillation. In fact, 
in the literature, the life time of the Rabi oscillation is of- 
ten taken as a measure of the coupling strength to the en- 
vironment. In this paper we point out a non-perturbative 
effect of the environment which can not be encapsulated 
in the level broadening picture. 

This effect results in a quantum phase transition as 
a function of the ratio of two coupling strength, 7/ A, 
when the cavity mode is in resonance with the two-level 
system. The numerator, 7, is the interaction strength 
between the two-level system and the cavity mode. The 
denominator, A, is same quantity between the two-level 
system and the environment. The phase at large 7/ A 
exhibits damped Rabi oscillation, while that at small 
7/ A does not. One can tune through this phase transi- 
tion by changing the degree of excitation of the cavity 
mode. This can be achieved by, e.g., varying the Rabi 
driving strength. The universality class of this phase 
transition is that of the inverse-square Ising model in one 
dimension. It is triggered by the proliferation of up <-» 
down flips of the two-level-cavity system. The purpose 
of this paper is to make predictions for the manifesta- 
tion of this phase transition in an absorption experiment. 

The model In the absence of the environment we use 
the Jaynes-Cummings modelQ to describe the coupling 
between a two-level system and a single cavity electro- 
magnetic mode, 



Here a z = ±1 denotes the upper/lower states of the two- 
level system, and a + (a) is the creation (annihilation) op- 
erator of the frequency ljo cavity photon. In Eq. (fT]) 
a± = a x ± itjy, and 2t$ is the energy gap of the two level 
system. 

The Hilbert space upon which Eq. (fT]) acts is the di- 
rect product of the Hilbert spaces of the two-level system 
and the cavity photons. A state in this Hilbert space is 
labelled by \s,n) where s = ±1 denotes the states of the 
two-level system and n counts the number of cavity pho- 
ton. It's easy to see that the matrix corresponds to H 
decouples into independent 2x2 blocks. Each block is 
spanned by the states | + , n — 1) and |— , n). If we use a 
pseudospin variable t z to label these two states so that 
\ Tz = +1) = 1 + 1, n-1) and \r z = -1) = |-l,n), Eq. © 
becomes 



H = A t z + -yr x , 



(2) 



where Ao = to — ^o/2. In obtaining the above equation 
we have dropped a constant energy term (n— l/2)o;o and 
defined 7 = 2g^/n. The quantity 7 measures the coupling 
between the two-level system and the cavity photon. It 
can be varied by changing g or n. Experimentally reach- 
ing large 7 by making n large can be achieved by coupling 
the two-level system to a strong (hence classical) driving 
field - the Rabi driving field. 

The environment is modeled by a bath of harmonic os- 
cillators in the spirit of Caldeira and Leggett In this 
paper we consider a specific type of environment whose 
coupling with the two-level-cavity system is given by 



H = A q t z + jt x + J2 9»T z {b a + bl) + ^2uj a blb a . (3) 



Here a labels the environment oscillators and u> a is their 
frequency. As shown by Caldeira and Leggett [l(| the 
information about the environment is summarized in the 
following spectral function J(uj) = ~^2 a g1 t 5{bJ — u> a ). In 
the rest of the paper we shall consider an "Ohmic" envi- 
ronment i.e., 



Hq = to<J z + <?(c+£i + <7-oJ) + ujoa'a. 



(1) 



J(co) — A u) e 



/u 



(4) 



2 



where 





(b) 
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FIG. 1: (a) The Rabi oscillation P(t) as a function of t/^/y) 
calculated using perturbation theory. In producing this graph 
we took 7 = 0.05to,A = 0.01, wo = 2t ,uj c = 2t . (b) P(t) as 
a function of tA (A is defined in Eq. (|19[) ) in the instanton 
metal phase. In producing this figure we took the Q = A in 
Eq. CE2I)- 



where A is a dimensionless coupling constant and u> c is a 
cut-off frequency. 

In perturbative picture the type of coupling in Eq. 
will give rise to dephasing. In the jargon of nuclear 
magnetic resonance the dephasing lift time T 2 measures 
the strength A in Eq. The model in Eq. ^ can be 
realized by, e.g., a Rabi driven "Cooper pair box" in tun- 
neling contact with a metal. The type of environmental 
effect described by Eq. §5§ can be simulated by the 
charge noise caused by quantum tunnelling of electrons 
between a trapping center in the tunnel junction and 
the metal. 
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Perturbative treatment Given that at t = the 
two-level system is in, say, the upper state, the Rabi os- 
cillation will be manifested as the quantum oscillation of 
the two- level-cavity system between r z = ±1. It can be 
seen by considering the following quantity 



P(t) = \(+l\e 



-iHt] 



l ) |2_| ( _ 1 | e -«ft| + 1) |2_ (5) 



For simplicity let us consider the case of resonant driv- 
ing, i.e., loq = 2t . In the absence of the coupling to 
the environment (g a — 0) it is simple to show that 
P(t) = cos(2jt), the Rabi frequency is 71-/7. To the low- 
est order, the retarded self-energy due to the coupling to 
the environment is given by 

E s , s (u) = {s'\V[u-H +ief 1 V\s) (6) 

where s,s' = ±1, V = J2 a 9aT z (b a + b^), and H = 
l T x + Y] a u a b' a b a . Standard calculation gives 



y 1 2 V/2M AM, 



(7) 



S 

f 2 (uj) = s IM^ + s l) - i7 rJ(u + sj)8(lu + s 7 )](8) 

s 

In the above equations 9{lo) is the step function, and 



A(w) = -A 



lu c + LOC 



dt — 
t 



Using the above result and the relation 



( S '\e- lHt \s) 



~2m 



u) + ie — Hq — 



(9) 



8), (10) 



we have calculated P(t) for a typical set of parameters, 
and the result is shown in Fig. 1(a). Indeed the Rabi 
oscillation is damped. 

Non-perturbative treatment We can integrate 
out the environmental oscillators to obtain the following 
quantum partition function for the two-level system [12j 



E 



dT 2n / dT 2n -i ■ 



dn 7 2 " e ~ s l Q ^ 
(11) 

Here 2n counts the number of t z — > — t z flips of the two- 
level-cavity system occurring between imaginary time 
t = and r = (3. The quantity Q(ti) is an alternat- 
ing sequence of ±1 (+1 for down — * up flip and —1 for 
up — > down flip) at ti, ..,T2 n . Setting Ti — ► 1 the action 
S[Q(Ti)} in Eq. ([TTJ) is given by 

s = J2Q i G( nj )Q j + A J2(Qi - - n), (12) 



where 
G(r) 



du> 



J(w) 



exp(-«|r|) = Aln(|r|/r ) (13) 



with To = l/w c - I n the limit of tr a << 1 and A tq << 1 
perturbative renormalization group calculation gives, 



dX 
dA 



dl 



A (l-7 2 ). 



dl \ 2 11 



(14) 



here I = ln(r/ro) is the running cutoff scale. For ujq = 2to 
these equations predict a phase transition between an 
instanton (i.e.,r 2 flipping event) insulator phase and an 
instanton metal phase. In the former phase 7 is renor- 
malized to zero at large time scale, while in the latter it 
diverges. 
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Alternatively one can map Eqs. (|llll2|) onto the classi- 
cal partition function of the inverse square ferromagnetic 
Ising model in one dimension 

? H = - K T,T^T2^ (15) 

where a is the lattice constant. By computing the Boltz- 
mann weight associated with domain wall configurations 
it can be shown that the lattice constant a in the Ising 
model plays the role of To and 

1 = e - ClK , \ = c 2 K, and A = 2h/a, (16) 

where c\ and ci are order unity constants. The advantage 
of mapping onto the inverse square Ising model is that 
for h = extensive Monte Carlo simulation has been 
performed on such model [l3| . This will allow us to access 
the physics in the tr® ~ 1 regime. 

Now let us briefly review the simulation results of 
Eq. (p~5|) at h = 0. As a function of K there are two 
phases. The weak K phase is a paramagnet and the 
strong K phase is a ferromagnet. Due to the long-range 
nature of the spin-spin interaction the phase transition 
occurs at finite K despite of one-dimensionality. In the 
paramagnetic phase the spin-spin correlation decays 

exponentially (oicr,,) ~ f | g . ° x .| J e~^ Xi ~ x ^^, and the 
correlation length diverges as £/a ~ e b/VK c -K ^ ^ e m 
the Kosterlitz-Thouless transition. The ferromagnetic 
phase has non-zero (cr,). At the ferromagnetic — » 
paramagnetic transition (<Tj) undergoes a discontinuous 
jump to zero. The value of (at) at the critical point 
is universal and is equal to 1/2. In the ferromagnetic 
phase the connected part of the spin-spin correla- 
tion function decays algebraically with distance, i.e., 
T(xi - Xj) = {diCJj) - (cr^icrj) ~ \ Xi - x -\n - Tne deca y 
exponent 77 is expected to be equal to 4y/ (K — K c )/K c . 
At the critical point T{x) ~ l/ln|-|. In the presence of 
a non-zero magnetic field, one expects (<7j) 7^ for all K 
values. The phase transition will not be present because 
the global <7j — > — <7j symmetry is explicitly broken. In 

this case one expects T(x) ~ f]fj) e~' x '/£. 

Predictions for optical absorption experiment 

In this section we use the results discussed in the last 
section to predict the outcome of optical absorption ex- 
periment. (Note that if strong 7 is achieved by Rabi 
driving, this would amount to a pump and probe exper- 
iment). Imagine coupling the two-level system to the 
following AC probe 

H pr obe = g c cos(ut)T z . (17) 

By applying the Fermi golden rule, it is simple to show 
that the absorption spectrum is given by 

A(u>) = g 2 c J2 \(a\r,\0)\ 3 S(w -E a + E ). (18) 



In the above |0) is the ground state and \a) are the ex- 
act eigenstates of the combined two- level- environment 
system. It can be shown simply that the zero- 
temperature imaginary time spin-spin correlation func- 
tion is the Laplace transform of A(u>): (t z (t)t z (0)) = 
/ °° duje~^ T A{yj) I g 2 c . Equivalently the absorption spec- 
trum is the inverse Laplace transform of (t z (t)t z (0)). 
Thus the knowledge of the behavior of the spin-spin cor- 
relation function in the last section enable us to make 
explicit predictions about the absorption spectrum. 

First we consider the case of resonance, i.e., ojq = 2io- 
This corresponds to the abscissa, 7/A, of Fig. 2. Ac- 
cording to Eq. (fl"6)) this is equal to e~ clK jc^K of the 
Ising model. Thus large/small 7/A corresponds to a 
small/large K. As a result we expect the system to be in 
the instanton metal phase for large 7/A and in the instan- 
ton insulator phase in the opposite limit. After inverse 
Laplace transforming the spin-spin correlation function 
we obtain the following prediction for the optical absorp- 
tion spectrum at low frequencies 

A(uj) = k 1 0(u>- A)(w- A)e- uj/n (19) 

in the instanton metal phase, and 

A(lj) = $ (7/A) 6(w) + k 2 uJ TI ~ 1 e- UJ / Q . (20) 

in the instanton insulator phase. In Eqs. p9[) and 
([20]) . fci,2 are constants, £1 is the absorption bandwidth, 
3>(7/A) = is the magnetization of the Ising model. 
Eq. (fl9|) predicts the presence of an absorption gap A 
in the instanton metal phase. This gap is related to the 
correlation length in the Ising model by A ~ ^ | . As the 
Ising correlation length diverges at the critical point A 
closes according to 

A^r - 1 exphV\/(7/A)-(7/A) c ]. (21) 

At criticality a zero-frequency delta function in the ab- 
sorption spectrum abruptly develops and the finite fre- 
quency absorption has a 1/u dependence. Throughout 
the instanton insulator phase, Eq. (|2T)|) predicts the pres- 
ence of a zero- frequency delta function peak and a gapless 
absorption spectrum. 

The case of off-resonance (Ao ^ 0) corresponds to the 
presence of a magnetic field in the Ising model. Translat- 
ing the behavior of Ising spin-spin correlation into A(ui), 
we obtain 

A(uj) = $ (7/A, A ) %) + k 3 8(uj - A){u - A)"- 1 e- w / n 

(22) 

Similar to the instanton insulator phase, a zero- frequency 
delta function is present. On the other hand, the non- 
zero frequency absorption shows a gap similar to the in- 
stanton metal phase. By fixing 7/A at a low value and 
tuning loq one will see the closing of the absorption gap 
when loq reaches the resonance value. 
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FIG. 2: A schematic phase diagram. The thick line in the 
abscissa is the line of critical points discussed in the test. 
All regions except the thin line in the abscissa are in the 
instanton insulator phase. A scan along the line with an arrow 
will exhibit the closing of the absorption gap when cross the 
abscissa. 



One might think that reaching the critical point on 
the abscissa of Fig. 2 requires the fine tuning of two pa- 
rameters. This is in fact not true, because the entire 
A = 0, < 7/A < (7/A) c interval corresponds to a line 
of fixed points. This is similar to the low temperature 
critical line of the Kosterlitz-Thouless phase transition. 
To reach this critical line one simply sets 7/A at a low 
value and adjusts the detuning parameter Ao- In terms 
of optical absorption, the critical line is reached when 
the finite frequency absorption gap closes. This enables 
one to tune to resonance (i.e., the abscissa of Fig. 2) eas- 
ily. Once this is achieved, one can tune across (7/A) c by 
adjust the Rabi driving strength. 

To study the Rabi oscillation in the instanton metal 
phase we need to compute the P(t) in Eq. ([5]) in the 
presence of an environment. To do it non-perturbatively 
we connect it to the correlation functions of the Ising 
model. Let |0) be the ground state of Eq. ([3]). The initial 
observation of t z — + 1 projects the state of the system to 
( 2* ) 1^)- Using the time evolution operator associated 
with Eq. to evolve this state and computing the t z 
expectation value at a later time t we obtain 



P(t) 



1 



{(0|r 2 (t)|0) + (0\T z T z (t)r z \0) 

(23) 



2(0|l + r z |0) 

(0\r z r z (t)\g) + (0\T z (t)T z \g)}, 



where T z (t) = e r z e~ . It is important to note that 
in the above equation t is the real (not imaginary) time. 
At Ao = and in the absence of spontaneous symmetry 
breaking (0|t*(*)|0) = (0|t z t z (£)t 2 |0) = 0. As a result 

nOC 

P(t) = (0\T z T z (t)\g)+(0\T z (t)T z \g) = / cos(wt)A(w)dw. 



(24) 

Inserting Eq. (fT9|) into Eq. (|24|). we calculated the Rabi 
oscillation in the instanton metal phase, the result is pre- 
sented in Fig. 1(b). One can see a damped Rabi oscilla- 
tions qualitative similar to the perturbative result shown 



in Fig. 1(a). In the instanton insulator phase, the strong 
coupling with the environment locks the the two-level- 
cavity system, and the Rabi oscillations is complectly 
quenched. 

In summary, we have shown that although the en- 
vironment is harmful to quantum coherence, it con- 
tains interesting physics in its own right. In particu- 
lar, we have demonstrated that when the environment is 
Ohmic and only causes dephasing, there is a quantum 
phase transition between a state which exhibits damped 
Rabi oscillation to a state which does not. In reality 
the dephasing-only requirement of our model has been 
achieved partially already in experiment. For example, 
for the "Cooper pair box" a T\ (life time due to the up 
«-> down flipping of the two-level system) which is an 
order of magnitude longer than T 2 (the dephasing life 
time) has been achieved 14|. In the langauge of Ising 
model, Ti/tq imposes a finite size cutoff to the critical 
behavior discussed earlier. Hopefully future technologi- 
cal progresses will enable one to probe the regime where 
T\ >> T2, hence the novel ID analog of the Kosterlitz- 
Thouless critical behavior. 
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